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We study Mott insulators of fermionic alkaline earth atoms, described by Heisenberg spin models 
with enhanced SU(iV) symmetry. In dramatic contrast to SU(2) magnetism, more than two spins 
are required to form a singlet. On the square lattice, the classical ground state is highly degenerate 
and magnetic order is thus unlikely. In a large- N limit, we find a chiral spin liquid ground state with 
topological order and Abelian fractional statistics. We discuss its experimental detection. Chiral 
spin liquids with non- Abelian anyons may also be realizable with alkaline earth atoms. 



An exciting thread in the study of strongly interact- 
ing cold atomic gases is the notion that such systems 
can be used as quantum simulators of strongly correlated 
materials Simple model systems can be engineered 
with a high degree of control, and studied as analogs of 
solid state materials. On the other hand, in some cold 
atom systems the simplest realizations of strong correla- 
tion physics may have no solid state analog. This raises 
the exciting prospect of systems and phenomena that are 
thus far unanticipated. 

Recently, it has been argued that fermionic alkaline 
earth atoms (AEA) in optical lattice potentials can real- 
ize a variety of model correlated systems, many of which 
lack solid state analogs and are relatively unexplored the- 
oretically Q • Fermionic AEA have nuclear spins as large 
as 1 = 9/2 for 87 Sr; due to lack of hyperfme coupling 
with the electronic ground state ( So), the nuclear spin is 
essentially decoupled from the electronic degrees of free- 
dom. This decoupling, also present in the lowest elec- 
tronic excited state ( 3 Po), implies that the s-wave scat- 
tering length is independent of nuclear spin, and leads 
to an enlargement of the spin rotation symmetry from 
SU(2) to SU(Af), where N = 21 + 1$,% This observa- 
tion, together with recent progress in and prospects for 
manipulating AEA opens the door to experimental 
studies of SU(iV) magnetism. We shall see here that the 
enlarged symmetry has striking physical consequences. 

In this Letter, we study the simplest antiferromagnctic 
square-lattice SU(iV) Heisenberg model that can be re- 
alized with AEA in the electronic ground state. We find 
that, as in some geometrically frustrated systems, for 
N > 3 magnetic order is underconstrained and there 
is a large degeneracy of classical ground states. Here, 
the degeneracy arises not from geometrical frustration 
but from the structure of the SU(iV) exchange interac- 
tion, and is present on any lattice for large enough N. 
This result indicates that magnetic order is unlikely, so 
we focus instead on non-magnetic ground states, which 
are controllably accessed in a large- N limit, where we 
find the ground state is the long-sought chiral spin liquid 
(CSL) [i, [SHU- The CSL spontaneously breaks time- 



reversal (T) and parity (V) symmetries, and is closely 
related to fractional quantum Hall liquids, sharing their 
remarkable topological properties Q . 

Specifically, we consider the large-?/ (insulating) limit 
of a Hubbard model with m < N atoms per site. N < 10 
can be realized with 87 Sr by populating a subset of the 
nuclear spin levels 2\. For m = 1, the spin at each site 
transforms in the fundamental representation of SU(iV), 
and N sites are needed to form a singlet, a crucial differ- 
ence from SU(2) magnetism. While m = 1 best avoids 
three-body losses, we also consider m = N/k for integer 
k > 2; in this case k sites are needed to form a singlet. 
Such models, which may be realizable for m not too large, 
allow us to consider a solvable large- N limit, where N is 
taken large with k fixed. This is a large- N generalization 
of the model with m = 1 and N = k, as the number of 
sites needed to form a singlet is preserved. 

It is convenient to define the model in terms of fl a 
(a = 1 . . . ,N), which creates a fermion on the square 
lattice site r. The Hamiltonian is 

H = J ]T S^{r)S^{r'), S af3 (r) = f r J rf) , (1) 

(rr') 

where the sum is over nearest-neighbor bonds, and J 
is the exchange energy. We have a local constraint, 
frafra = m - Study of correction terms arising away 
from the large-?/ limit will be deferred to future work. 

Most studies of SU(JV) magnetism have focused on 
models where two sites can be combined to form a sin- 
glet. The most-studied cases are the k = 2 model defined 
above [l(|, and models defined by placing conjugate rep- 
resentations on the two sublattices of a bipartite lattice 
Spin-3/2 alkali fermionic atoms exhibit an enlarged 
SO (5) symmetry, where also two sites can be combined to 
form a singlet [l2j. Finally, we note that the models we 
discuss here have been solved exactly in one dimension for 
rn = 1 13]. In two dimensions, the N = 4, m = 1 model 
has been studied in the context of orbitally-degenerate 
Mott insulators, although there the SU(4) symmetry re- 
quires substantial fine-tuning [3]. On the cubic lattice, 
plaquette states (see Fig. [TJ;) have been studied using a 
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quantum plaquette model [15| . 

Semiclassical limit. The semiclassical limit considered 
here is a generalization of the large-S 1 limit of SU(2) mag- 
netism. We consider a generalized model where the spin 
at each site transforms in the SU(iV) irreducible represen- 
tation labeled by the Young tableau with one row and n c 
columns [H| . This representation is the symmetric com- 
bination of n c fundamental representations, and in the 
SU(2) case is a spin-5 spin (S — n c /2). 

We can define this model in terms of fermion opera- 
tors /J aa , where a = 1, . . . , n c is a "color" index. On 
every site we place n c fermions, and antisymmetrize over 
their color indices. Defining S a p(r) = J2 a fraafrpa' the 
Hamiltonian is identical in form to Eq. ([I]). We define 
the coherent state \z) = (z a fai) ■ ■ ■ {z a fl n J\0), which 
is parametrized by the A^-component complex spinor z 
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z only changes \z) by a 



phase, the overall phase of z is unphysical and coher- 
ent states are labeled by points in the manifold CP — 1 , 
which has dimension 2(N — 1). In the limit n c — > oo, 
the state Y[ r \ z r) is an eigenstate, and the energy is 
E = Jn 2 c J2(rr')\4^\ 2 + 0(n c ). 

The energy is minimized for z\z , — on nearest- 
neighbor bonds. For N > 2, we immediately see a signif- 
icant difference from SU(2) magnetism: knowing z r does 
not uniquely determine the neighboring z r < that mini- 
mizes the energy. This leads to an extensive degeneracy 
of classical ground states. To see this, we estimate the di- 
mension D of the ground state manifold [16j | . Letting N s 
be the number of lattice sites, the total dimension of all 
the CP^ 1 spins is 2N S {N-1). On every bond, z\z r , = 
provides two constraints, for a total of AN S constraints. 
Treating the constraints as independent leads to a lower 
bound: D > 2N S (N - 3). For N = 3, where this bound 
is not helpful, it can be shown by explicit construction of 
ground states that D oc N s . 

Such extensive degeneracy is a hallmark of geometri- 
cally frustrated systems, where underconstraint emerges 
from the inability to simultaneously satisfy a set of com- 
peting interactions. A crucial physical consequence is 
a strong, even complete, suppression of magnetic order 
[ill. The semiclassical limit is biased towards magnetic 
order, and since it is suppressed even there, we expect 
that the present models lack magnetic order altogether 
for n c = 1, the case of interest for AEA Mott insulators. 

Large-N limit. Returning to the model Eq. |T]), mag- 
netically disordered ground states can be be controllably 
studied in the limit N — > oo, where m — N/k, J = J '/N , 
and k and J are held fixed. This limit was studied for 
k = 2 in 1 1(1. w here the ground state is a valence-bond 
solid (VBS) [17(. Here, we find the k — 3, 4 ground states 
break lattice symmetry and are analogous to the VBS 
(Fig. QJ. For 5 < k < 10 we present evidence that the 
ground state is the CSL, and also discuss low-lying com- 
peting states. We conjecture that the CSL is the ground 
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FIG. 1: Large- N dimer and plaquette ground states for k — 2 
(a), k = 3 (b) and k — 4 (c). Xrr' has constant magnitude on 
the dark bonds and is zero on the others. For k = 3 (k — 4), 
the phase of \ rT t is chosen so that the flux through each 
plaquette is 7r (zero). The patterns shown are not necessarily 
those selected by 1/N corrections. 



state for all k > 5. 

The mathematical structure of the large- N limit is the 
same as for the k = 2 case already studied. The prob- 
lem reduces to finding the ground state of the mean- 
field Hamiltonian Hmft = Hmft + J2 r ~ /™/™). 
where Hmft = (N/J)J2{rr') \Xrr>\ 2 + U K and H K = 
J2( rr ')(Xrr'fi a f r > a + H.c). This is required to satisfy 
the self-consistency conditions 
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(#«/ra> (a). rn=UtJ ra ){b). (2) 



The field Xrr' arises from decoupling the exchange inter- 
action on each bond, and fj, r arises from a Lagrange mul- 
tiplier field implementing the constraint of m fermions 
per site. Without loss of generality, we assume J2 r fi r = 
0. A set of (Xr-r'jMr) satisfying Eq. ([2]) is a mean-field 
saddle point. The saddle-point energy -Emft is an ex- 
ternum with respect to variations of the fields, but not 
necessarily the global minimum. The task at hand is to 
find the lowest energy saddle point as a function of k. 

For k = 2, Rokhsar established a lower bound on 
-Emft, and showed that, on any lattice where a dimer 
covering is possible, any dimer state such as that shown 
in Fig. QJi saturates the bound [ItJ • The leading correc- 
tions in the 1/N expansion then select an ordered VBS 
configuration from among the various dimer states 11 1. 
It is straightforward to extend Rokhsar's bound to gen- 
eral k. First, for a given saddle point, using Eq. ((2Jd), 
Emft = (Hmft) > Emft, the ground state energy of 
Hmft- A lower bound on Emft is easily obtained fol- 
lowing Ref. [ItJ. For k = 2, one divides the spectrum 
of Hk in half; in general, one divides the spectrum into 
occupied and unoccupied levels. On the square lattice, 
one finds .Emft > ~[(k - l)NJN s ]/2k 2 . 

A stricter lower bound can be established for bipartite 
lattices, where the spectrum of TLk is symmetric about 
zero energy. We divide the spectrum into the sets £ 
(occupied levels), U (image of C under e — * — e), and M 
(remaining levels). An analysis similar to that of Ref. 17 1 
shows that, on the square lattice, 



Emft > -NJN s /4k. 
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FIG. 2: Lowest-energy competing state for fc = 6 (the fc = 5 
state has a similar pattern). The lattice is covered by stripes, 
of which one is shown. The shading of bonds represents 
\Xrr' I) interpolating between the maximum \xrr' I (black), and 
|Xrr'| = (white). Some regions enclose 7r-flux, as indicated. 



fc 


CSL 


ICSL 


LC 


5 


-0.043080 


-0.043070 


-0.042987 


6 


-0.033069 


-0.03299 


-0.032961 


7 


-0.026130 


-0.02597 


-0.025730 


8 


-0.021138 


-0.02102 


-0.020897 



TABLE I: Energies of CSL and competing states, in units of 
NJNs = N 2 JN S , for 5 < fc < 8. ICSL is the lowest-energy 
inhomogeneous CSL that was found. LC is the lowest com- 
peting state that cannot be interpreted as an inhomogeneous 
CSL. Note that the energy difference between CSL and LC is 
larger for fc = 7, 8 than for fc = 5, 6. 

For fc > 2 this bound is stricter than that above, and is 
saturated if and only if J2 r ^rn r — 0, C (U) contains only 
the constant energy — e (+e), and M. contains only zero 
energy levels. Here, h r is the average fermion number 
calculated in the ground state of Hk ■ 

For fc — 3, 4 the bound is saturated by the plaque- 
tte states (Fig. [IJ. In each case 1/N corrections are ex- 
pected to select an ordered configuration, and the large- 
N ground state for k = 3,4 is analogous to the k = 2 
VBS. It is impossible to saturate the bound for k > 8: the 
large- N ground state energy of the two-site problem is 
-NJ/k 2 , which gives the bound £mft > -2NJN s /k 2 , 
stricter than Eq. ([3]) for k > 8. Even for 5 < k < 8, the 
conditions needed to saturate Eq. ([3]) are very restrictive 
and we conjecture they cannot be satisfied. 

Below, we present evidence that the CSL is the large- 
N ground state for 5 < k < 10. The CSL saddle point 
has iJ. r = and |Xrr'| = X- with the phase of \rr' such 
that the flux through each square plaquette is 2ir/k. This 
results in a fermion band structure with k bands, where 
only the lowest band is filled, and for fc > 3 it is separated 
from the others by a gap. This mean-field state is a lattice 
integer quantum Hall state: there is a quantized Hall 
conductance of N , where the (fictitious) fermion charge 
and Planck's constant have been set to unity j 1 81 ] - 

To determine the lowest-energy saddle point for 5 < 
k < 10, in addition to explicit construction of saddle 
points, we implemented a numerical self-consistent min- 
imization (SCM) algorithm. The algorithm begins with 
a random choice of Xrr', and self-consistently iterates 
equation while choosing pL r at each step to satisfy 



equation d2ji>) . We allowed Xrr' and fi r to vary within 
a given unit cell embedded within a larger system (with 
periodic boundary conditions). It can be proven that 
SCM converges to a local minimum of the energy. For 
5 < k < 8, we studied all rectangular unit cells with fc 2 
or fewer sites, excluding cells of unit width. For each cell, 
we ran the SCM procedure on at least 30 (in some cases 
more than 500) different sets of random initial conditions. 
The CSL was the lowest energy state found (TableQ]). For 
fc — 9, 10, less extensive application of SCM also found 
no states below the CSL in energy. 

We also find locally stable competing states, some only 
slightly higher in energy than the CSL (Table []}. The 
competition between CSL and these states will need to 
be resolved by going beyond the large- iV limit, and, ul- 
timately, by experiments. The lowest such states found 
can be viewed as inhomogeneous versions of the CSL. For 
fc = 5, a 2 x 2 ordering pattern is superimposed on an 
average 27r/5 flux per plaquette. For fc = 6, 7, 8, the CSL 
divides into domains. As long as the CSL remains stable 
to inhomogeneity (e.g. the domain wall energy is posi- 
tive), these states will not be ground states. Therefore, 
we also searched for the lowest competing states that can- 
not be viewed as inhomogeneous CSLs. For fc = 5, 6, we 
find stripe states (Fig. [2]) that break various lattice sym- 
metries but preserve T. For fc = 7, 8, we find a distinct 
CSL with 2n/2k flux per plaquette. 

Properties of CSL. The CSL is characterized by both 
its broken symmetries and topological order. T and V 
breaking is signaled by a nonzero spin chirality (C123) 7^ 
0. Here, C123 = i(Pi23 ~ -P321) is the spin chirality of 
lattice sites 1,2,3, and P123 the operator that cyclically 
permutes the spin quantum numbers on those sites @ . 

Understanding topological order requires going beyond 
the mean-field description. It is important at this stage to 
note that /J Q does not create an atom. Instead it creates 
a spinon, which carries the spin but not the conserved 
atom number. The most important fluctuations about 
the saddle point are in the phase of \rr f ~ {Xrr / )c ia " rr ' •, 
where the spatial component of a fluctuating U(l) 

gauge field coupled to the spinons. The time-component 
of the gauge field arises from the fluctuations of fj, r . 
The gapped spinons can be integrated out to obtain a 
Chern-Simons (CS) effective action for the gauge field, 
S c f{ — (N/4n) J dtd 2 r a^e^xd^ax, where the coefficient 
is determined by the mean-field Hall conductance. The 
CS term is responsible for the topological properties of 
the CSL 0. It converts spinon excitations into anyons 
with statistical angle it + ir/N. Moreover, its presence 
implies the spinons are deconfined and propagate freely, 
and the CSL thus exhibits quantum number fractional- 
ization. For a system with an edge, there are gapless 
chiral edge modes. Finally, the ground state degeneracy 
on a surface of genus g is 2 A'' 9 , where the factor of 2 arises 
from the spontaneous T-breaking. 

Experimental detection. The distinct features of the 
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states discussed here can be split into two categories: 
straightforward ones associated with the spin gap and 
broken symmetry, and, for the CSL, more subtle prop- 
erties having to do with the presence of topological or- 
der. A number of well-developed experimental techniques 
can be employed to detect the features of the first type. 
Radio-frequency spectroscopy can be used to see the pres- 
ence of the gap [L9J . One can measure spin-spin correla- 
tion via noise correlations to see the absence of order [2(| • 
The VBS-analog states for k = 3, 4 could be detected by 
adiabatically merging groups of sites into a single site, 
followed by application of bandmapping techniques [2l[ . 
To detect T-breaking, one can superimpose a second sys- 
tem of fermions ( 3 Pq alkaline earths [2]) or bosons (al- 
kali atoms) , that couples to the CSL atoms via spin-spin 
interaction. By symmetry, this coupling will induce an 
effective orbital magnetic field for the second system. For 
bosons, this field will induce vortices, and for fermions it 
will lead to detectable changes of the energy spectrum. 

Topological order detection is more challenging. Due 
to the chiral edge modes, a disturbance of the spins near 
the system edge will propagate around the edge with a 
well-defined velocity and direction, which could poten- 
tially be detected. Returning to the Hubbard model, the 
CSL will exhibit spin-charge separation. Letting c£, cre- 
ate an atom, there will be a particle carrying the atom 
number but not the spin, created by b' = ct/ a . The 
6-particle can be thought of as a bound state of an atom 
and a spinon, and has fractional statistics with angle 
tt/N. These expectations can be formalized using a slave- 
rotor treatment of the Hubbard model [22j • Because the 
spinon does not couple directly to a scalar potential, low- 
ering the optical potential at a lattice site can localize a 
6-particle, the fractional statistics of which could poten- 
tially be probed by techniques proposed in the context 
of quantum Hall- like states in cold atomic systems [23| . 

To observe these characteristic properties, the temper- 
ature should be at most on the order of the gap A to the 
lowest-energy quasiparticle excitations. Using the large- 
N limit and boldly setting N = k (i.e. one atom per 
site) , in the CSL we find A ~ J for both excitations of the 
gauge field, and particle-hole excitations of the spinons. 
The harmonic trapping potential determines the spatial 
extent of m — 1 Mott insulator; provided this is larger 
than the characteristic scales of the CSL, its signatures 
can be observed. Using the large- N limit, these length 
scales are estimated to be at most a few lattice constants. 

Finally, we note that the n c = 2 model, discussed in the 
context of the semiclassical limit, can be realized using 
one ground state atom and one 3 Pq atom on each site, 
depending on the sign of the Kondo exchange 0. We 
have shown, and will present in detail elsewhere, that 
this model can support a CSL with a fluctuating U(2) 
gauge field, with a SU(2) CS term at level N. This CSL 
supports non-Abelian anyons, and is a candidate system 
for universal topological quantum computation [24 1. 



In summary, we studied square lattice Mott insulators 
that can be realized by fermionic AEA in the electronic 
ground state. We showed that magnetic order is unlikely, 
found the CSL ground state in a large- TV limit, and dis- 
cussed its experimental detection. 

We thank Paul Julienne for useful discussions, and ac- 
knowledge support from NIST and NSF (A.M.R.) and 
NSF grant no. DMR-0449521 (V.G.). 
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